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Abstract 



Oh 



We define the quadratic algebra su(2) Q which is a one-parameter deformation of the Lie 
algebra su(2) extended by a parity operator. The odd-dimensional representations of su(2) 
i-Jh I (with representation label j, a positive integer) can be extended to representations of su(2) a . 

We investigate a model of the finite one-dimensional harmonic oscillator based upon this al- 
\ gebra su(2) a . It turns out that in this model the spectrum of the position and momentum 

operator can be computed explicitly, and that the corresponding (discrete) wavefunctions can 
be determined in terms of Hahn polynomials. The operation mapping position wavefunctions 
into momentum wavefunctions is studied, and this so-called discrete Hahn-Fourier transform is 
computed explicitly. The matrix of this discrete Hahn-Fourier transform has many interesting 
properties, similar to those of the traditional discrete Fourier transform. 

m 

1 Introduction 

t-H ! 

' In standard theory of quantum mechanics, position and momentum operators are (essentially) self- 
adjoint operators in some infinite-dimensional Hilbert space, satisfying the canonical commutation 
relations. Quantum mechanics in finite dimensions has attracted much attention in recent years pQ. 
In a finite-dimensional Hilbert space, the canonical commutation relations no longer hold. Despite 
this, finite-dimensional quantum mechanics has been useful in areas such as quantum computing 
and quantum optics PIE]. The defining relations for a quantum mechanical oscillator in a finite- 
dimensional Hilbert space are not unique [I], and several type of models have been proposed. Our 
interest comes mainly from models related to some Lie algebra (or a deformation thereof). The finite 
oscillator model that has been studied most extensively is based on the Lie algebra su(2) or so (3) in 
the one-dimensional case [SHZ] • In the case of a two-dimensional oscillator this has been generalized 
by the same authors to so (4). Such finite oscillator models are of particular importance in optical 
image processing [5], and more generally in models where only a finite number of eigenmodes can 
exist. For example signal analysis dealing with a finite number of discrete sensors or data points led 
to physical models realizing a one-dimensional finite oscillator [8rll0| . In a previous paper , a new 
model for the finite one-dimensional harmonic oscillator was proposed based on the algebra u(2) Q , a 
one-parameter deformation of the Lie algebra u(2). This u(2) a model offers an alternative position 
and momentum spectrum compared to the su(2) model. Furthermore, the position wavefunctions 
have simple expressions in terms of Hahn polynomials, with interesting properties related to those 
of a parabose oscillator 
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For a one-dimensional finite oscillator, one considers three (essentially self-adjoint) operators: a 
position operator q, its corresponding momentum operator p and a (pseudo-) Hamiltonian H which 
is the generator of time evolution. These operators should satisfy the Hamiltonian-Lie equations 
(or the compatibility of Hamilton's equations with the Heisenberg equations): 

[H,q]=-ip, [H,p]=iq, (1) 

in units with mass and frequency both equal to 1, and h = 1. Furthermore, one requires [6]: 

• all operators q, p, H belong to some (Lie) algebra (or superalgebra) A; 

• the spectrum of H in (unitary) representations of A is equidistant. 

The case with A = su(2) (or its enveloping algebra) has been treated considerably in a number 
of papers [SHI]- The relevant representations are the common su(2) representations labelled by 
an integer or half- integer j. In such a representation, the Hamiltonian is taken as H = J$ + 
j + \, where Jo = J z is the diagonal su(2) operator. Thus the spectrum of H is n + | (n = 
0, 1, . . . , 2j). With q = \{J+ + J_) = J x and p = §(J+ — J-) = —Jy, the relations (P) are satisfied. 
Clearly, q and p have a finite spectrum given by {—j, — j + 1, . . . , +j} [6]. More important, the 
position wavefunctions have been constructed, and are given by Krawtchouk functions (normalized 
symmetric Krawtchouk polynomials). These discrete wavefunctions have interesting properties, 
and their shape is reminiscent of those of the canonical oscillator [6]. This is explained by the fact 
that under the limit j — > oo the discrete wavefunctions coincide with the canonical wavefunctions 
in terms of Hermite polynomials [6|ll2j. 

In [11], the case with A = u(2) Q was investigated (a > —1). For this one-parameter deformation 
of u(2), only the representations labelled by a half-integer j survive as representations of u(2) Q (so 
only the even-dimensional representations). This led to an alternative model of the finite oscillator, 
with the spectrum of H again given by n + \ (n = 0, 1, . . . ,2j). For the operators q and p, the 
spectrum is again finite and equidistant in steps of one unit, except that there is a gap of size 2a + 2 
in the middle; explicitly, it is given by 
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- a- j - - j + -,■ ■ ■ ,-a - l;a + l,a + 2,. . . ,a + j + -. (2) 

An interesting result in [11] was that the position wavefunctions could be constructed and they 
turned out to be normalized Hahn (or dual Hahn) polynomials. An investigation of these discrete 
wavefunctions gave rise to remarkable plots, suggesting in fact a relation with the parabose oscilla- 
tor. Indeed it was shown that under the limit j — > oo the discrete wavefunctions coincide with the 
parabose wavefunctions in terms of Laguerre polynomials |llj . 

Despite the novel results presented in [11] . it remains somehow intriguing that the odd-dimen- 
sional u(2) representations can not be deformed as u(2) a representations. In the present paper, we 
present the solution to this problem. It turns out that one should consider a different one-parameter 
deformation of su(2), involving only the three su(2) operators and a parity operator (and no central 
element). This new deformation, denoted by su(2) Q , is defined and it is shown that in this case the 
odd-dimensional su(2) representations, labelled by an integer value j, can be deformed as su(2) Q 
representations. Using su(2) a as a model for the finite oscillator, yields the same spectrum of H 
given by n + \ (n = 0, 1, . . . , 2j). The finite spectrum of the position and momentum operators q 
and p is however quite different, and given by 

0, ±y/k(2a + k + l), (k=l,...,j). 

The position (and momentum) wavefunctions are constructed, and turn out to be again Hahn 
polynomials (in this case with parameters (a, a) or (a + 1, a + 1)). These wavefunctions are once 



2 



more discrete versions of the parabose wavefunctions, as a limit computation shows. A fascinating 
question in this context is the relation between the position wavefunctions and the momentum 
wavefunctions. In the canonical case, these functions are related through the Fourier transform. 
Here, we show that these wavefunctions are related through a discrete version of the Fourier trans- 
form, which we refer to as the discrete Hahn-Fourier transform. For readers primarily interested 
in special functions, the computation of this discrete Hahn-Fourier transform is probably the most 
attractive part. It involves a special case of a bilinear generating function (the Poisson kernel) for 
dual Hahn polynomials. 

The contents of the remaining sections is as follows: in section 2 the deformed algebra su(2) Q 
and its representations are constructed. In section 3 we use su(2) Q as a model for the finite 
oscillator, and determine in particular the spectrum of the position and momentum operators, and 
their eigenvectors. The structure of these eigenvectors is studied in section 4, yielding position 
and momentum wavefunctions. In section 5 we determine the operation that transforms position 
wavefunction into momentum wavefunctions, the so-called discrete Hahn-Fourier transform F, and 
its properties. The somewhat technical but interesting proof of the 4-^3(1) form of the matrix 
elements of F is given in section 6. Finally we give a summary of the results in section 7. 



2 The algebra su(2) a and its representations 

The Lie algebra su(2) [13l[I3] can be defined by its basis elements Jo, J+ , J- with commutators 
[Jo, J±] = ±J± and [J + , J_] = 2 Jo. The non-trivial unitary representations of su(2), corresponding 
to the star relations Jq = Jo, jj_ = J=p, are labelled [TBTH] by a positive integer or half- integer j and 
have dimension 2j + 1. The action on a standard basis vectors \ j, m) (with m = — j, — j + 1, . . . , + j) 
is given by 

Jo \j, m) =m \j, m) , J± \j, m) = =F m)(j ± m + 1) \j,m±l). 

The universal enveloping algebra of su(2) can be extended by a parity operator P with action 
P\j, m) = (— \y +m \ j,m). This means that P commutes with Jo, anticommutes with J+ and J_, 
and P 2 = 1. This extended algebra can be deformed by a parameter a, leading to the definition of 
su(2) a . 

Definition 1 Let a be a parameter. The algebra su(2) a is a unital algebra with basis elements Jo, 
J+, J_ and P subject to the following relations: 

• P is a parity operator satisfying P 2 = 1 and 

[P, J ] = PJ - J P = 0, {P, J±} = PJ± + J±P = 0. (3) 



The su(2) relations are deformed as follows: 



[J ,J±] = ±J±, (4) 
[J+,J-] = 2J + 2(2a + l)J P. (5) 



Note that this is different from the u(2) Q deformation introduced in [11] ; in particular it does 
not involve a central element. For the deformation u(2) a , the even-dimensional u(2) representations 
(j half-integer) could be deformed. The current su(2) a is "complementary" in the sense that now 
the odd-dimensional su(2) representations (j integer) can be deformed: 
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Proposition 2 Let j be an integer (i.e. 2j is even), and consider the space Wj with basis vectors 
\j,—j), | j, — j + 1), • • ., \j,j)- Assume that a > — 1. Then the following action turns Wj into an 
irreducible representation space o/su(2) Q . 



P\j,m) = (-iy +m \j,m), 
Jo\j> m ) = m \3i m ), 



J+\j, m) 
J- 1 j, rn) 



V (j - m)(j + to + 2a + 2) \j,m + 1), 
- to + 2a + l)(j + m + 1) |jf, m + 1) 

a/( j + m + 2a + 1) (j - m + 1) |j, m - 1) 
V 7 0' + - m + 2a + 2) |j,m - 1), 



if j + m is even; 

if j + to is odd, 

if j + m is odd; 

if j + to is ei>en. 



(6) 
(7) 

(8) 
(9) 



The proof is essentially by direct computation, the same as |114 Proposition 2]. Note that the 
representation given in this proposition is unitary under the star conditions = P, Jq = Jo, 
J_j_ = J=p. Also note that for a = — ^, the deformation is trivial (both for the algebra and the 
representations) . 

3 Using su(2) a as a model for the one-dimensional oscillator 

Quite similar as in the non-deformed case [5] or in the u(2) a deformed case [11], let us choose the 
position, momentum and Hamiltonian (representation dependent) operators as follows: 



i(J+ + J_), p=^(J+-J_), H = J +j + -. 



(10) 



These operators satisfy ([T]). In the representation space Wj, H\j,m) = (m + j + to), therefore 
the spectrum of H is linear and given by 

1 



n + 



(n = 0,l,...,2j). 



(11) 



From the actions (H])-©, one finds 



2q\j, to) = \J (j + to) (j — to + 2a + 2) | j, to 
if j + m is even, and 

2q\j, m) = \J {j — m + 1) (j + to + 2a + 1) | j, to 



1> + ^{j -m){j + m + 2a + 2) \ j,m + l), 
1) + + m + l)(i - m + 2a + 1) \ j, m + 1) 



if j + to is odd. The action of 2ip is similar. In the (ordered) basis {\j, — j), \ j, — j + 1), . . . , \j,j 
1)) li>i)} °f Wj, the operators 2q and 2ip take the matrix form 



2q 



2ip 
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M 










\ 




M 





Mi 
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M 2J . 


-i o 


1 


/ 
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with 

M = /V (fc + l)(2j + 2g-fc + l), if k is odd; 

k \ y / (k + 2a + 2)(2j - k), if k is even. 

For these matrices the eigenvalues and eigenvectors can be constructed in terms of Hahn polynomi- 
als. Hahn polynomials Q n (x; a, P, N) [T5l[16] of degree n (n = 0, 1, . . . , N) in the variable x, with 
parameters a > — 1 and P > — 1 are defined by [I~5|[l~6]: 



r> t a at\ t? f-n,n + a + P + l,-x \ 

Q n (x;a,P,N) = 3 F 2 \ a + 1 _ N ' 1 J' ( 15 ) 

in terms of the generalized hypergeometric series 3F2 of unit argument |17|,I18]. Their (discrete) 
orthogonality relation reads [T5|[l6]: 

JV 

} w(x; a, p, N)Qi (x; a, P, N)Q n (x; a, P, N) = h(n; a, ft, N) 8i n , (16) 



x=0 



where 



n a r\ fa + x\ (N + B — x\ 
w(x;a,P,N) = I ^ W R P _ ^ J (x = 0,1, . . . , N); 

w a AT\ {n + Ot + P + I) N +1 {P + l) n n\ 

h[n;a,P,N) — 



(2n + a + /3 + l)(a + l) n (JV - n + l) n JV! ' 

We have used here the common notation for Pochhammer symbols |17yi8| (a)k = a(a + 1) • • • (a + 
A; — 1) for /c = 1,2,... and (a)o = 1. Orthonormal Hahn functions Q are determined by: 



Q n (x; a, P, N) = = . (17) 

^h(n;a,P,N) 

Recall that dual Hahn polynomials have a similar expression. In fact, for x 6 {0, 1, . . . , JV} the 
right hand side of ([To]) is the dual Hahn polynomial R x (X(n); a, P, N) of degree x in the variable 
X(n) = n(n + a + P + 1); see [T5j[16] for their orthogonality relations. 

In terms of these, the eigenvalues and (orthonormal) eigenvectors of M q can be constructed: 

Proposition 3 Let M q = 2q be the tridiagonal (2j + 1) x (2j + I) -matrix (fl~2l) and let U = 

(Uki)o<k,i<2j be the (2j + 1) x (2j + 1) -matrix with matrix elements: 

U 2 i, j = (-l) i Qo{r,a,a,j), i G {0, 1, . . . , j}; £7 2mj = 0, i € {0, . . . ,j - 1}; (18) 
(-1V - 

U 2 i,j-k = U 2 i,j+k = —^-Qk{i;a,a,j), i € {0,1,..., j}, fc e {1, . . . ,i}; (19) 

(-1)* ~ 

U 2i+ ij-k = -U 2i+ ij + k = -j=-Qk-i(r,a + l,a + l,j - 1), i 6 {0, 1, . . . , j - 1}, 

ke{i,...,j}. (20) 

T/ien U is an orthogonal matrix: 

UU T = U T U = I, (21) 
the columns of U are the eigenvectors of M q , i.e. 

M q U = UD q , (22) 
where D q = diag(eo,ei, • • • ,^2j) is a diagonal matrix containing the eigenvalues of M q : 



ej^k = -2^k{2a + k + l), e d = 0, e j+k = 2y/k(2a + k + 1), (k = 1, . . . (23) 



5 



Proof. Using the orthogonality of the Hahn polynomials, and the explicit expressions ()18p - (|20p . 
a simple computation shows that (U T U) k i = 5 k i. Thus U T U = I, the identity matrix, and hence 
UU T = I holds as well. 

Now it remains to verify (I22D and (I23D . By the form of M q : 

2j 

{MiU). k = M u U ik = M^xUi-y, + MiU i+hk . (24) 
1=0 

We have to consider (|24p in six distinct cases, according to i even or odd, and to k belonging to 
{0,1, ... ,j — 1}, to {j + 1, j + 2, . . . , 2j} or k = j. Let us consider the case that i is odd and 
k € {0, 1, . . . , j — 1}. Then, relabelling the indices appropriately, (fl4"|) . (fT9j) and (f20j) yield: 

(M q U)2i+ij-k = M 2i U 2iJ _ k + M 2i+1 U 2 i+2, j -k 

= (-iyV2^(a + i + i)Q k (i; a, a,j) + (-l) i+1 V2^(i + i)(j + a - i)Q k (i + 1; a, a,j) 



At this point, the forward shift operator formula for Hahn polynomials |15|. (9.5.6)] can be applied 
and yields 



/„,„ TT x / -.4 /- (a + l)i+i(a + fe(fc + 2a + 1) ^ 

(M^) 2i+1 ,_, = (-l)^J ^-_ i ^)IMib;a;a, J -) (a + l)i « + ^ « + ^ ~ ^ 



-2v / Mfc + 2a + l)C/ 2i+lii _ fc = ej_ k U 2i+ld „ k = (UD q ) 



2i+lj-fc' 



For the other five cases, the computations are similar and for some of them the backward shift 
operator formula for Hahn polynomials [151 (9.5.8)] is applied. □ 
Note that in the case of u(2) Q [II] . the equation corresponding to (j24"|) was related to two new 
difference equations for Hahn polynomials [HUGO]. Here, the equation is just corresponding to 
known forward or backward shift operator formulas. The above proposition yields in particular the 
spectrum of the position operator: 

Proposition 4 The 2j + 1 eigenvalues q of the position operator q in the representation Wj are 
given by 

- y/j(2a + j + 1),-VU ~ 1)(2« + j), ■ ■ ■ , -V2cV+2; 0; V2c7+2, y/j(2a +7+1). (25) 
It will be appropriate to label these q-eigenvalues as q k , where k = —j, — j + 1, . . . , +j, so 



q± k = ±^k(2a + k + l), k = 0,l,...,j. 



One can compare the spectrum of q with that in related models, see Figure [IJ In the non-deformed 
case (a = — |, or the su(2) model), the spectrum is just equidistant as already mentioned in the 
introduction. In the case of u(2) Q (and j half-integer), the spectrum is given by ([2]): apart from an 
extra gap of size 2a + 2 in the middle, it is again equidistant. In the current case of su(2) a (and j 
integer), the spectrum given by Proposition 0] is not equidistant. 

Essentially the eigenvectors of (112p have components proportional to Hahn polynomials with 
parameters (a, a) when the component has even index and with parameters (a + l,a + 1) when 
the component has odd index. It is convenient to introduce a notation for these eigenvectors: the 
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orthonormal eigenvector of the position operator q in Wj for the eigenvalue q k , denoted by \j, q k ), 
is given in terms of the standard basis by 

3 

\j,Qk)= ^2 U 3+-m,j+k\j,m). (26) 

m=-j 

Let us now turn our attention to the momentum operator p. Up to signs, the matrix M p is the 
same as the matrix M q . The analysis of its eigenvalues and eigenvectors is thus very similar. We 
present just the final result here: 

Proposition 5 Let M p = 2ip be the tridiagonal (2j + 1) x (2j + 1) -matrix (|13p and let V = 

(Vki)o<k,K2j be the (2j + 1) x (2j + 1) -matrix with matrix elements 

V k ,i = -i k+l U kU (27) 
where U is the matrix determined by ()18p -([20 p . Then V is a unitary matrix, 

yyt = y]y = /. 

The columns ofV are the eigenvectors of M p , i.e. 

M'PV = VD p , (28) 
where D p is a diagonal matrix containing the eigenvalues e k of M p : 
D p = diag(e ,£i,. . . ,e 2 j), 



£j- k = -2iy/k{2a + k + \), ej = 0, e j+k = 2i v / fc(2a + k + 1) (k = 1, . . . , j). (29) 
Hence the 2j + 1 eigenvalues p of the momentum operator p in the representation Wj are given 

by 

- \J j(2a + j + 1),-VU ~ 1)(2« + j), • • • , -v^T2; 0; V2cV+2, ^/j{2a+j + 1), (30) 

in other words, the momentum operator p has the same spectrum as the position operator q. We 
shall denote these p-eigenvalues by p k , where k = —j, + so 

p± k = ±y / k(2a + k + l), k = 0,l,...,j. 

The normalized eigenvector of the momentum operator p in Wj for the eigenvalue p k , denoted by 
\j,Pk), is then given by 

3 

\j,Pk)= ^2 V 3+rn,j+k\j,m). (31) 
m=-j 

Finally, it is worth mentioning another property of the matrix elements of V, that follows from 
the explicit expressions (127p . (I18p - (l20p and the orthogonality properties of the Hahn polynomials: 



V T V 



( ••• -1 \ 
••• -1 



(32) 



\ -l ••• o o j 

Note also that from (I27D it follows that 

V = JU where J = diag(-i, 1, i, -1, . . .), (33) 

and the sequence (— i, 1, i, — 1) is repeated, ending with i or —i (since the matrices are odd- 
dimensional). One can also write J = — i diag(i°, i , i 2 , i 2 , . . . , i 2j ). 
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4 Position and momentum wavefunctions and their properties 

The position (resp. momentum) wavefunctions of the su(2) a finite oscillator are the overlaps be- 
tween the g-eigenvectors (resp. ^-eigenvectors) and the iJ-eigenvectors (or equivalently, the Jo- 
eigenvectors \j,m)). Let us denote them by $^; m (g) (resp. Vl/^ m (p) ), where m = + 
1, . . . , +j, and where q (resp. p) assumes one of the discrete values q k (resp. p k ) (k = —j, —j + 
1, . . . , +])■ Therefore, in the notation of the previous section: 

$j+ m (<?fc) = {j, m\j, q k ) = U j+m j +k , (34) 
^j+m(Pk) = (j, m \j,Pk) = V j+mJ+k . (35) 

Let us consider the explicit form of these wavefunctions, first for the position variable. For j + m 
even, j + m = 2n, and for positive g-values one has 

(a) (— l) n ~ 

$ 2n (Qk) = — 7=~Qk(n;a,a,j), n = 0, k = l,...,j, (36) 



V2 



or equivalently: 



( a ) (-l) n w(n;a,a,j) f-k, k + 2a + 1, -n 

®2n \Qk) = — 7^~\ Tn T 3 2 I , i • ;i 

V2 V h{k;a,a,j) V o> + l,-j 



(~l) n w(n;a,a,j) 
y/2 V h(k;a,a,j) 



Rn(ql;a,a,j) , g| = fc(2a + fc + 1), (37) 



where i? n (A&; a, a, j) is a dual Hahn polynomial [15^16] of degree n in the variable X(k) = q\ = 
k{2a + k + 1). In a similar way, one finds for j + m odd, j + m = 2n + 1, and for positive g-values, 
that 



(a) , _ (-1)" / w(n;a + l,a + l,j -1) f -k + l,k + 2a + 2,-n 

2n+AQk)- ^ V h (k-l;a + l,a + l,j-l) 3 * 2 \ a + 2, -j + 1 ' 



/ w(n;a + + 1, j - 1) , 2 , . ^ 

= V /.(fc-l;a + l,a + l, J -l) ^^- 2(a + 1);a + 1 ' a + 1 ' J - 1 )- (3§) 

For g = or negative values, the expressions are of course analogous, according to (fT8"j) - (f2"U]) . 

Before studying some special properties of these position wavefunctions, and determining the 
momentum wavefunctions, let us consider plots of these functions for some a-values. We choose 
a fixed value of j, say j = 30, and plot some of the wavefunctions &l?\q) for certain values of 
a. Since a = —\ is a special case (where su(2) a reduces to su(2)), there are three cases to be 
considered: — 1 < a < — |, a = —\ and a > — |. In Figure [2] we take a = — |, a = —0.7 and 

a = 2 respectively. We also plot in each case the ground state ^^\q), some low energy states 
^^(q) and $2*\q), and the highest energy state. The plots are similar as in the u(2) Q case (where 
j is half- integer) . The most obvious difference is that is part of the spectrum now. Another 
difference, but more difficult to see in the plots, is that the spectrum is not equidistant for a ^ 
For a = — ^, these plots coincide with the ones given in the su(2) model [5,6j. The wavefunctions 

are Krawtchouk functions. This follows also from the following transformations of 3F2 
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series to 2-^1 series when a = — |, according to [U (48)] 

f-k,k,-n \ U n(2n) w {-2n,-j-k 



3 r 2 



0) 

2i 



-fc + l,fc + l,-n A (-1)" L+l) 771 f-2n-l,-j-k \ 
3A -J + 1 2fe (i-i) ~2j 



(39) 



(40) 



2fc (V) 

For a/-^, the plots are comparable with the parabose wavefunctions [11] . One can indeed 
again study the behaviour of the discrete wavefunctions (q) when the representation parameter 
j tends to infinity. In this process, one should pass from a discrete position variable q to a continuous 
variable x. This can be done by putting q = j x l 2 x and then compute the limit of jV^nHq) for 
j — > 00. The actual computation is similar to the limit computation performed in [11] §4], so we 
shall not give any details. Note that, due to (|37|) . 



q 2 = (k + a+ l -f-{a+ l -)\ or * = -a - I + yjtf + (a + \f. (41) 

Using this last expression for k in the 3F2 expression of (|37|) . and replacing herein q by j 1 / 2 ^, the 
limit can be computed and yields: 



lim j^U^x) = (-l)\ r( u \xr^e^l-Lt\x\ (42) 

j^oo W T(a + n+l) 



in terms of Laguerre polynomials. Similarly, one finds from f|38|) : 



lim j 1 '**^ 1 ^) = (-l) n J r( , ^ sM^e^L^V). (43) 

j^oo y 1 [a + n + 2) 

The functions in the right hand side of (|42p - (|43p are known: they are the wavefunctions of the 
parabose oscillator [2TH23] with parameter a = a + 1 > 0, see the appendix of [TTJ. 

It remains here to consider the momentum wavefunctions. Due to the fact that the spectrum 
of q and p is the same, and due to the similarity of the matrix of eigenvectors V (compared to U), 
the expressions are analogous and we give only the final result here: 



= -^J h^a,a,j) Rn {ph *' *' J) ' Pl = + k + l) ' (44) 



^(p*) = y ^^1-1) ^ (ft - 2(a + 1); a + 1, a + l,j - l) . (45) 

5 The discrete Hahn-Fourier transform 

In canonical quantum mechanics, the momentum wavefunction (in L 2 (M)) is given by the Fourier 
transform of the position wavefunction (and vice versa): 



¥(p) = -L J e-^(q)dq. 
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In the current case, we are dealing with discrete wavefunctions, and we should look for an analogue 
of this. In terms of the notation of the previous section, let 



*(<fc) 



*i(«fc) 



*(p*) 



/ *o(p fe ) \ 

*l(p*) 



(k 



(46) 



So it is natural to define the discrete Fourier transform in this case as the matrix F = (Fik)-j<i,k<+j 
relating these two wavefunctions. In other words: 



*(p0 = Fm 

fc=-j 



(47) 



As this generalized discrete Fourier transform maps Hahn polynomials into Hahn polynomials, we 
shall refer to it as the discrete Hahn- Fourier transform. By (|34p - (|35p . the columns of V consist of 
the column vectors ^(pk) (k = and similarly for the matrix U. So (|47p actually means 

that V = UF, or: 

F = U T V. (48) 

Using the explicit matrix elements from U and V, (fl~8P - ([20j) and ([27]) . this leads to following form 
of the matrix elements of F: 



F j-k,jTi = F j+k,j±i = -^^2(-^) n Qk(n;a,a,j)Qi(n;a,a,j) 

n=0 

± 2 J2(-l) n Qk-i(n) a + 1, a + 1, j - (n; a + 1, a + 1, j - 1) 

n=0 

• j 

F 3Tk,j = F j,jTk = — -^^2(-^) n Qk(n;a,a,j)Q (n;a,a,j), k = l 

* n=0 

-i^(-l) n Qg(n;a,a,j). 



(49) 
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k,l = l. 



(50) 



(51) 



n=0 



In the following section, we shall determine explicit expressions for the above matrix elements. 
But before that, we can already summarize some properties of the discrete Hahn-Fourier transform 
matrix F. 

Proposition 6 The (2j+l) x (2^+1) -matrix F is symmetric, F T = F , and unitary, F^ F = FF^ = 
I. Furthermore, it satisfies F 4 = I, so its eigenvalues are ±1, dbi. A set of orthonormal eigenvectors 
of F is given by the rows of U , determined in Proposition [3 The multiplicity of the eigenvalues 
depends on the parity of j. When j = 2n is even, then the multiplicity of —i, 1, i, —1 is n + 1, n, n, n 
respectively. When j = In + 1 is odd, then the multiplicity of —i, —1 is n + 1, n + 1, n + 1, n 
respectively. 



Proof. The symmetry of F is easily seen from the expressions (I49p - (l5ip . The unitarity of F follows 
from (|48p . the orthogonality of the real matrix U and the unitarity of V. Again using (|48p and 
the orthogonality of U, one finds F 2 = F T F = V T UU T V = V T V. But the explicit form of V T V 
is known, see (f32j) . Since (V T V) 2 = I, the result F A = I follows. So the eigenvalues can only be 
±1, ±i. Using again (051) and (I33J), one sees that F = U T V = U T JU, or 



FU T = U T J. 
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j-n 



In other words, the columns of U T (or the rows of U) form a set of orthonormal eigenvectors of 
F, and the eigenvalues of F are found in the diagonal matrix J . From the explicit form of J , the 
statement of the multiplicities follows. □ 

Note that these properties are similar to those of the common discrete Fourier transform [9|I24| . 
For the case a = — 5, the matrix F was already studied in [9]. In that special case, the position and 
momentum wavefunctions (|46p are in terms of Krawtchouk polynomials, and the corresponding 
discrete Fourier transform can be referred to as the discrete Krawtchouk-Fourier transform. The 
matrix elements of F are in that special case again Krawtchouk functions [9]. 

One of our main results is the explicit computation of the elements of F for general a. Note 
from P9|) -([5T ]) that all these matrix elements are of the form Yln=o(~ ^) n Qk(n] <^> &,j)Qi( n ] a, a,j), 
with < k, I < j. So apart from a factor {hlk^a^a^^h^^a^,^))' 1 / 2 , this expression is equal to: 

3 

>.,j)Ql(n;a,a,j) (52) 
Q k (n;a,a,j)Qi(n;a,a,j) (53) 
R n (X(k); a, a, j)R n (X(l); a, a,j). (54) 

J-n J 

In the form (|52p . the right hand side is just like the orthogonality relation of Hahn polynomials, 
except for the extra factor (— l) n . In the form (f54l) . one can see that S(k,l,a,j) is a special case 
of the Poisson kernel (or a bilinear generating function) for dual Hahn polynomials: 

E* n (° I ") ( P y_ ~ U ) R n( x (k); a, p,j)R n (\(iy, «, P,j). 

However, as far as we know a closed form expression for this Poisson kernel is not known. In fact, the 
best one can do is re-express the product R n (X(k); a, (3, j)R n (X(l); a, (3,j) into a 8-P7 hypergeometric 
series (following the approach of [25, §8.3], where the g-analogue is given). Here we will show that 
this Poisson kernel does have a simple expression when /3 = a and t = — 1. 

We shall prove the following result, yielding the explicit expression for the discrete Hahn- Fourier 
transform matrix F: 

Theorem 7 The special Poisson kernel for dual Hahn polynomials S(k,l,a,j), where k and I are 
integers with < k, I < j, satisfies 

S(k,l,aJ) = (-l) k+l+ iS(kJ,a,j), (55) 



S(k,l,a,j) = 




) n w(n; a, 


n=0 






3 

= £(- 


-i; 


]U /a + n 


71=0 




\ n 






3 

= E(- 




]n fa + n 


n=0 




V n 
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so it is whenever k + 1 + j is odd. For the other cases, it is given by: 



S(2K, 2L, a, 2J) /(i^M'^^W Wtl), 

x M~ K ' K Vltrh + l +l A < 56 > 
\ a + J + 1, a + 1, — J J 

2 2J (\)j- K (\)j_ L {a + + J + 2) K (a + J + % 

(2J)!J(i) J 

x 4F3 f-^ + « + i-^ + « + § ;1 ^ (57) 



5(2^ + 1,2^ + 1,0,27) 



a + J + 2,a + l,-J + l 

S(2K, 2L + l,a,2J + l) 



2 2J+1 {\)j- K+ i{\)j-L{a + 1)j+i(« + J + 2)x(a + J + 2) L 



(2J + l)!(i) J+1 

\ a + J + 2, a + 1, — J J 

and finally S(2K + 1, 2L, a, 2 J + 1) is given by replacing K and L in the right hand side of (jSHJ). 

First of all, note that S(k, I, a,j) = for k + I + j odd implies that in expression ff^9|) only one of 
the two parts survive (either the real part or else the imaginary part). Together with Theorem [7] 
this implies that each element Fn- is, up to a factor, equal to a terminating Saalschutzian 4^3(1) 
series. In other words, up to a factor each F^ is also a Racah polynomial [15]. The unitarity of 
the matrix F does not lead to any new relations for Racah polynomials, but just follows from their 
discrete orthogonality relations. 

Before discussing the proof of Theorem [TJ let us examine what happens to this discrete Halm- 
Fourier transform in the limit when j — > 00. Suppose that j = 2J is even, and let us first consider 
the limit of the imaginary matrix elements of F. When k and I are even (j = 2J, k = 2K and 
I = 2L), then according to (fi9l) : 

F j+k,j±i = ~(h(k; a, a, j)h(l; a, a, j))~ 1/2 S(2K, 2L,a,2J), 

with S(2K, 2L, a, 2J) given by (|56p . Just as for the limit computation of the wavefunctions, see (|4ip . 
it is necessary to make the replacements 



k - \/q 2 k + (a + i) 2 - (a + i), 1= ^ pf + (a + ^) 2 - (a + i), 

and put qk = j l / 2 x and p\ = i x l 2 p in the above expression. After doing this, the limit of the 4-F3 
series appearing in S(2K, 2L, a, 2J) is fairly easy to determine by termwise computation: 

v 77 ( ~ 2 ' 2 a ^~ 2 ' ~~ 2 ' 2 a ~^ 2 r 1 tp f ~ 1 2 2 
J^oo V a + 1 + 1, ck +- 1, — § / V a + 1 4 

The last expression is proportional to J a (xp) , where J a is the Bessel function of the first kind [26 1 . 
The limit of the factors in front of the 4-F3 series are elementary but a bit more tedious to compute, 
and we shall not give the details of this computation. Similarly, one has to determine the limit of 
the real matrix elements of F. Adding both contributions, one finds: 

lim j 1/2 F j+kJ± i = - (\xp\ 1/2 J a (\xp\) + ixp\xp\~ 1/2 J a+1 (\xp\)) . 
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The function in the right hand side is known: it is the kernel of the generalized Fourier transform, 
studied by Mukunda et al in the context of the parabose oscillator |21j . So one can conclude that 
our discrete Hahn-Fourier transform tends to the generalized Fourier transform of [21] for large 
values of j. 



6 Computation of the discrete Hahn-Fourier transform matrix el- 
ements 

The purpose of this section is the proof of Theorem [71 In fact, we shall show (|55p and (|56p ; the 
proof of ([57]) and (158]) is analogous. 

Let us first collect some known transformation formulas for hypergeometric series. The first is 
Thomae's transformation of a terminating Saalschiitzian 4-^3(1) series |18|. (2.4.1.7)], [171 (7.2.1)]: 

F /a, b, c, -AT \ = (/ - c)jy(e + f-a-b) N p f e-a,e-b,c,-N \ 
4 V e J,9 J (/);v(e + / - a-b- c) N 4 3 \e,e + f - a - b,e + g - a - b' J 

where e + / + g = 1+a+b+c— N. For the Hahn polynomials appearing in (152j) . the following 
transformation formula holds [27, p. 186]: 

Q k (x; a, a, N) = (-l) k Q k (N - x; a, a, N). (60) 

In the same paper [27\ (3. 18)- (3. 19)], one can find the following 4.F3 expressions for Hahn polyno- 
mials with parameter (3 = a: 

„ , ^ ( —k,k + a + 7y,—x,x — N \ » » 

^ a ' a ^=^(a;i,V,(l-JV)/2 ;1 J' (61) 

Q 2k+1 (x;a,a,N) = ^^ 4F J -^ + a + \ -x,x-N \ 

^ +n ' ' ' ; N \a + l,(l - N)/2, (2-N)/2 J v ; 

Next, we need an equation that goes back to a classical identity for 6j-coefficients. For 6j- 
coefficients of su(2), the following holds [281 (9.8.4)]: 



(63) 



It turns out that we need the corresponding identity for 6j-coefficients of positive discrete series 
representations of su(l, 1), which reads: 

\ > U k 3 ,k,k 23 U k 2 ,k,k 2 3 ~ U k 3 ,k,k!3 ■ ^ ' 

323 

Herein, U^^ 2 k ^ 2 is the standard notation for the Racah coefficient of su(l, 1) [29]. In (|64l) . k\, hi 
and ks are su(l,l) representation labels (i.e. they are positive real numbers), and in a common 
notation [29j 

fci2 = h + k 2 + 312, k = ki2 + k + 3 + /, k 13 = h + k 3 + j 13 , k 23 = k 2 + h + j 23 , 

with j'12, f, j'13 and j'23 all nonnegative integers. The summation in the right hand side of (|64[) 
runs over all j'23 with < j'23 < j' + ji2- Using the explicit expression in terms of a 4^3(1) series 
for the Racah coefficients in (I64p . given by [291 (4.41)], and making the following replacements: 

(/, ji2, J13) — >(p,q,r), (y,y,y) — > (a,b,c), j 23 — > n, 
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one arrives at the following identity between terminating Saalschiitzian 4^3(1) series: 

(p + q\ (b + c + 2n - 1) (b + c) n (a + b + c + p + g - l) n 
£^ V n / (6 + c + ra-l) (6 + c + p + q) n (l - a -p - q) n 

f—n,n + b + c—l,—q,q + a + b—l \ ( —n,n + b + c— l,—r,r + + c— 1 
x 4-T3 I ; 1 I 4-T3 I ; 1 

V 0, a + 6 + c + p + g— 1,— p — g / V c,a + b + c + p + q—l,—p — q 

= , 1 xp- r (fr + c Wg (Q)r(a) g 4j?3 / -9,9 + + 6- l,-r,r + o + c-l > A 
(o)p +g (6) g (c) r 43 Vo,a + 6 + c + p + g-l,-p-g'y 

This is not yet in the form needed for our proof of Theorem [71 But after performing Thomae's 
transformation ([59]) on the 2nd and 3rd 4^3(1) series appearing in ([65]) . and replacing r by p + q — r, 
we have the following: 

Lemma 8 Let p, q and r be nonnegative integers, and a, b and c arbitrary parameters, then 

« J (b + c + n-1) {c) n {b + c + p + q) n 

f—n,n + b + c—l,—q,q + a + b—l 
X 4-T3 I ; 1 

\ b, a + b + c + p + q— 1,— p — q 

( —n, n + b + c — 1, — r, r + 1 — a — c — 2p — 2g 

x 4-F3 , ; l 

V b,l-a-p-q,-p-q 



-i-iy 



{b + c) p +g(l - c - p - q) r 



(c)p(a + 6 + c + p + g- l) g (l -a-p-q) 



r 



P /-g,g + o + fe-l,-r,r + l-o-c-2p-2g \ 
x , n ; 1 . 66 

V b, 1-c-p-q, -p-q ) 

Note that the identities (I65p - (]66p show at first sight some similarity with the expansion formulas 
constructed in [30] (or some g-analogues in [31]); they turn out to be quite different however. 

We now come to the final part of this section: 
Proof of Theorem First of all, from (|52p . note the symmetry S(l,k,a,j) = S(k,l,a,j). 
Starting from (153|) . using (160p for Qk and Qi, and then reversing the order of summation (i.e. 
replace n by n — j), one finds (f55|) . This implies that S(k, I, a,j) = whenever k + l + j is odd. So 
we need to simplify the expression only when k + I + j is even. We shall do this explicitly in the 
case that k, I and j are even: k = 2K, I = 2L and j = 2J, with K, L and J nonnegative integers 
(the other three cases are similar). So we need to compute 

S(2K, 2L, a, 2J) = ^(-lf \ + J (° \f_ ~ n ) Q2 K (n; a, a, 2 J)Q 2L (n; a, a, 2J). (67) 

n=0 V n / V n J 

Using ([6T]) . this can be written as S(2K, 2L, a, 2 J) = X^n=o ^n(K, L), where 

/—if, K + a+ \,—n,n- 
\ a + l,-J,-J+i 

- 2J \ 

;lj. (68) 

Making an appropriate replacement of the summation variable, it is easy to see that 

2J j 
Y,Tn(K,L) = J2T n (K,L), 

n=J n=0 





14 



hence we can split the total sum YlnLoT n (Ki L) in two equal parts; taking care of the overlapping 
middle term, there comes 

S(2K, 2L, a, 2 J) = 2 [ J] T n (K, L) + \tj{K, L) \ . (69) 

\n=0 / 

In order to perform the summation in the right hand side of (j69[) . one can make the following 
substitution in ([66|) : 

p = J — K, q = K, r = L, a = -, b = a + t, c = —2J — a, (70) 

and then take the limit t — > 1 (in fact, for all terms in (|66p one can immediately take t = 1, except 
for the term with n = p + q = J the limit process is actually necessary). Consider, after the 
substitution (|70p . the nth term in the left hand side of (|66p. with 0<n<p + q = J: the coefficient 
in front of the two 4^3(1) 's becomes: 

(2 J)! , N „ fa + n\ fa + 2 J - n 
-(-1) 



(a + 1)2J V n ) V 2 J - n 

Consider similarly the term with n = p + q = J\ here the coefficient becomes, when t — )■ 1: 

1 (2J)! . ,jfa + J\fa + J 



2(a + l) 2 j V J J\ J 

Furthermore, it is easy to see that under this substitution (|70p and t — > 1 the two 4^3(1) 's in the 
left hand side of ([66]) become the 4^3(1) expressions of (f68|) . Thus, up to the coefficient ? 
the right hand side of (|69p can be summed according to (|66p . and one finds, after simplifications: 



S(2K, 2L,a,2J) 



2 2J (i)j_^(i) j_ L (a + l)j(a + J+l) K (a + J+ l) L 



(2J)!(|)j 

/ —if, if + a + I 

4^3 

proving ([56]) . □ 



, . -K,K + a + \,-L,L + a + ± 
x 4^3 I T ~ z ;l 

a + J + 1, a + 1, — J 



7 Summary 

The most popular finite oscillator model, especially for applications in quantum optics, is based on 
the Lie algebra su(2) or so (3) [SHI]- The dimension of the model depends on the su(2) representation 
label j, which is an integer or half-integer positive number. Its mathematical properties have been 
studied well, in particular the properties of the wavefunctions given by Krawtchouk functions. 
These wavefunctions have interesting plots, and in the limit j — > 00 these discrete wavefunctions 
tend to the continuous canonical oscillator wavefunctions. The discrete Fourier transform turning 
position wavefunctions into momentum wavefunctions has also been investigated in this case [9]. 

A first type of deformation of this su(2) model was offered by its q-deformation. The su q (2) 
model for the finite oscillator was investigated in [32,33]. The main properties are: the position 
operator has a discrete anharmonic spectrum, and the wavefunctions are given in terms of dual 
q-Krawtchouk polynomials [33]. 

We have now considered two different type of deformations of the su(2) model, closely related 
to each other. The first deformation u(2) Q was given in [11] and allows a deformation of the even- 
dimensional representations only {j half- integer). The second deformation su(2) a was the topic 
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of this paper, and allows a deformation of the odd-dimensional representations only (j integer). 
Both cases have a deformation parameter a > — 1 and for a = — | they reduce to the nondeformed 
su(2) model. The (discrete) spectrum of the position operator can be constructed explicitly in the 
deformed algebras. In the case of u(2) a it is equidistant except for an extra shift in the middle of the 
spectrum. In the case of su(2) a it is no longer equidistant but has the simple form ±^k(k + 2a + 1) 
(k = 0,1, . . . For k sufficiently large, this behaves like ±(k + a + ^), so it tends to an equidistant 
distribution for large k. 

The position and momentum wavefunctions have been constructed explicitly for the deformed 
models. They are given in terms of Hahn polynomials. For u(2) Q , the even wavefunctions are nor- 
malized Hahn polynomials with parameters (a, a + 1), and the odd wavefunctions with parameters 
(a + l,a). For su(2) a , the even wavefunctions are normalized Hahn polynomials with parameters 
(a, a), and the odd wavefunctions with parameters (a + 1, a + 1). The plots of these discrete wave- 
functions have nice properties, and in both cases they tend to the parabose wavefunctions when j 
is large. For a — > — ^ they tend to the Krawtchouk wavefunctions of the su(2) model; and of course 
for a — y — 7j and j — > oo they tend to the canonical oscillator wavefunctions in terms of Hermite 
polynomials. 

An interesting extra object studied in this paper is the operation that transforms position 
wavefunctions into momentum wavefunctions, i.e. the discrete analogue of the Fourier transform. 
For the case of su(2) a , this discrete Hahn- Fourier transform has been constructed explicitly, and is 
determined by a matrix F. This matrix shares many classical properties with the common discrete 
Fourier transform matrix. In fact, it has the extra interesting feature that there is a natural basis 
of eigenvectors of F. The main computational result of the paper is the proof that the matrix 
elements of F have a simple form in terms of terminating Saalschiitzian 4^3(1) series, i.e. in terms 
of Racah polynomials. Note that for u(2) Q this Hahn- Fourier matrix was not determined in |11| . 
but knowing the results for su(2) a it should be a routine computation to do this. 
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Figure 1: Plots of a typical spectrum of the position operator, (a) in the case of the su(2) model, 
(b) in the case of the u(2) Q model, and (c) in the case of the su(2) Q model. 
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0.4 

n=60, a=-m 




Figure 2: Plots of the discrete wavefunctions <3?n (q) in the representation with j = 30. The four 
top figures are for a = —1/2, the middle figures for a = —0.7, and the four bottom figures for 
a = 2. In each case, we plot the wavefunctions for n = 0, 1, 2, 60. 
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